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Introduction 

In previous work we have developed a novel approach to visualizing second 
order symmetric 2-D tensor fields based on degenerate point analysis. At 
degenerate points the eigenvalues are either zero or equal to each other, 
and the hyperstreamlines about these points give rise to trisector or wedge 
points. These singularities and their connecting hyperstreamlines determine 
the topology of the tensor field. In this study we are developing new methods 
for analyzing and displaying 3-D tensor fields. This problem is considerably 
more difficult than the 2-D one, as the richness of the data set is much 
larger. Here we report on our progress and a novel method to find , analyze 
and display 3-D degenerate points. First we discuss the theory, then an 
application involving a 3-D tensor field, the Boussinesq problem with two 
forces. 

Definitions 

Second-Order Tensor field 

Let C ( X , Y) be the set of all the linear transformations of the vector space X 
into the vector space Y, and let E be an open subset of R n . A second-order 
tensor field T (x) defined across E is a mapping: T : E — ► £ (R m , R m ) that 
associates to each vector x of E a linear transformation of R m into itself. If 
R m is referenced by a Cartesian coordinate system, T (x) can be represented 
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by m 2 Cartesian components T tJ (x) ,i,j = that transform according 

to 

T'ij = PipPjqTpq ( 1 ) 

p,q = \ 

under an orthonormal transformation ft = {ftij} of the coordinate axes. 1 

2-D and 3-D tensor fields can be represented by 2 x 2 and 3x3 matrices 
respectively. In Cartesian coordinates the matrices representing symmetric 
tensor fields take the following forms: 


• 2-D symmetric tensor fields 


T (x) = T (x,y) = 


( T n (x,y) 
\ T u (x,y) 


Tu(x,y) \ 
T 2 i{x,y) ) 


(2) 


• 3-D symmetric tensor fields 


T(x) = T(x,y,z) = 


( T u (x,y,z) 
Tu ( x,y,z ) 
Tu(x,y,z) 


T n (x,y,z ) T \3 (x, y, z) \ 
T 2 2 (x,y,z) T 2 3 (x,y, z) 
T 23 {x,y,z) T 33 (x,y,z) ) 


(3) 


Degenerate Points 

A degenerate point of a tensor field T : E — ♦ £(R. m ,R m ), where E is an 
open subset of R m , is a point x 0 £ E where at least two of the m eigenvalues 
of T are equal to each other. 

In the case of 2-D tensor fields, x 0 is a degenerate point iff A! (x 0 ) = 
A 2 (^o)- F° r 3-D tensor fields, various types of degenerate points exist, 
corresponding to the conditions Aj (2:0) = A 2 (^o)> A 2 (xo) = A3 (xo), or 
Ai (x 0 ) = A 2 (x 0 ) = A 3 (fo)- 


Locating Degenerate Points 

• 3-D tensor 

The first thing to note about a 3-D symmetric tensor field (Eq. (3)) is 
that it has 6 independent variables, three of which are on its diagonal. As a 
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result, various types of degenerate points may exist. These types correspond 
to the following conditions: 


Aj (x 0 ) = A 2 (x 0 ) > A 3 ( x 0 ) 
Ai (x 0 ) > A 2 (x 0 ) = A 3 (x 0 ) 
Ai (xq) = A 2 (xq)' = A 3 (xq) 


( 4 ) 

( 5 ) 

( 6 ) 


A degenerate point that corresponds to the condition described by Eq. ( 6 ) 
implies that the tensor at that location always appears in its diagonalized 
form: 

/ A 0 0 \ 


T(x 0 ) = 


0 

Vo 


A 

0 


0 

A ) 


( 7 ) 


The current method for locating this type of degenerate point involves solving 
5 simultaneous equations. 


' Tu(x 0 ) ~ T 22 (x 0 ) = 0 
T 22 (x o) - r 33 (x 0 ) = 0 

< r 12 (x 0 ) = o (8) 

Ti 3 (xo) = 0 
. T 23 (x 0 ) = 0 

Since we are working in a 3-D space, grid cells are cubes (Fig. (1)). Let /, 
g, p, q and r be the values of T\\ — T 22 , T 22 — T 33 , T 12 , T 13 and T 23 at the 

vertex (i,j,k) and f max , J min, 9max, Smim Pmaxy Pminj 9moi» r man r min 

be the maximum and minimum values at these eight vertices respectively. 



(i, j+1, k+1) 


(i+l, j+1, k+l) 




(i.j, k) 


(i+l. j. k) 



Figure 1: A 3-D grid cell. 


Only when all of the following conditions: 

' fmin < 0 

/max > 0 
9min — 0 
Qmax ^ 0 
Pmtn ^ 0 
Pmax ^ 0 

(/min — 0 
( Qmax ^ 0 
I'min — 0 
I’mai ^ 0 


\fmin | 

| + 

| fmax 

l> 

0 

Iffmin | 

+ 


l> 

0 

\pmin | 

| + 

| Pmax 

l> 

0 


+ 

| Qmax 

l> 

0 

|^m»n | 

; + 

\r max 

l> 

0 


( 9 ) 


are met, will there possibly be a degenerate point inside the cell. Then 
the next step is to interpolate f,g,p,q,r to approximate the tensor within 
the cell. In order to do the linear interpolation , we need to do the first 
order approximation for all of the f,g,p,q,r which gives us 5 corresponding 
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variables s,t,u,v,w. At this moment the physical meaning associated with 
the two extra variables s and t is unclear, where as the other 3 parameters 
represent the local position variables u,v,w. 

We can readily see that merely checking whether or not all 15 conditions 
are satisfied for each cube is a lot of work. Also, we always have to solve 
5 nonlinear equations simultaneously. Furthermore, this method cannot be 
applied to degenerate points that correspond to the conditions described by 
Eqs. (4,5). In these cases the tensor may look like: 


( A, 


T(f 0 ) = 


V 


0 

0 


( 10 ) 


or 


Ai 0 0 \ 

T(f 0 )=| 0 A 2 0 (11) 

0 0 A 2 / 

However, these tensors will be diagonal matrices only when they are in their 
eigenvector space. More generally, in any other spaces, they will take on the 

form: t u (x) r 12 (x) t 13 (x)\ 

T (x) = T u (x) T 22 (x) T 23 (x) (12) 

\T 13 (x) T 23 (x) T 33 (x\) 

A continuous tensor field varies at different locations, but we only have 
sampled data at grid points. Therefore we won’t be able to know the eigen- 
vector space of a field point off the grid, we usually choose a space by inter- 
polating grid points in its neighborhood, but the tensor will appear just as 
a regular tensor. It is very difficult to find a simple rule to check its double 
degeneracy that will work in any random space. 

Because of these difficulties, we went back to the original tensor and 
looked at it in a different way. 

Since a second-order tensor representation is a matrix, we can get its 
eigenvalues by solving its characteristic equation. 


• 3-D case 

A 3-D symmetric tensor field can be expressed as a 3 x 3 matrix(Eq. (3)) 
and its characteristic equation is: 

A (A) = —A 3 + aA 2 + bX + c 


(13) 




Figure 2: A typical three-root polynomial of order 3. 


where 


a. — T\\ + T 22 + T33 

b = T? 2 + Tl 3 + T \ 3 - r„r 22 - t„t 3 3 - T 22 T 33 (14) 

c = T u T 2 2 T 33 + 2T n T l3 T 23 — T u Tl — T 22 TI ^ — T 33 T 12 


and it’s derivative is: 


dA (A) 
d \ 


— — 3A 2 ■)■ 2aA -f- b 


(15) 


The characteristic equation A (A) is a polynomial of order three and as 
such it can be generally plotted as in Fig. (2). The number of roots repre- 
sents the number of different eigenvalues at each location: three-roots means 
non-degeneracy, two-roots means double degeneracy, one-root means triple 
degeneracy. 

The various possible cases of degeneracy are determined by the solutions 
to Eq.(15). 

. a ± \fd 

A„„ = — 5 — (16) 


where 


d = a 2 + 36 

The three possible degenerate cases are as follows: 
1. A; = Aj = A 2 > A 3 


(17) 
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Figure 3: Two-root polynomial of order 3. 



Figure 4: Two- root polynomial of order 3. 
A(Xi) = 0 


2. Ai > A 2 — A3 — A,, 

3 . A 1 = A 2 = A3 = A 1 = A„ 

Also, 


A (A„) = 0 


A (A, = A„) = 0 

d ( A,' = A,,' ) = a? T 36 = 0 


(18) 


( 19 ) 

( 20 ) 
( 21 ) 
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Figure 5: One-root polynomial of order 3. 


After the solutions to the roots of the derivative equation being substi- 
tuted into the characteristic equation, we get three spatial functions: 



2a 3 + 9 ab + d 3 ^ 2 


B\ (x,y,z) = A{ A,) = 

[■ c 

27 

(22) 


2 a 3 + 9ab-d 3 ' 2 


B 2 (x, y, z) = A(A«) = 

h c 

27 

(23) 

5 3 (x,y,z) = d(A, 

= A,j) = a 2 + 3b 

(24) 


These are the conditions for both double degeneracy and triple degener- 
acy. Our problem has been reduced to finding points where the functions 
B\ (x,y, z) = 0, B 2 (x, y, z) = 0, or B 3 ( x , y, z) = 0. 


Locating the Singular Points 

In 3-D space, B\ (x, y,z) = 0 is a maximum for B x (x, y, z), B 2 (x, y, z) = 0 is 
a minimum for B 2 (x, y, z) and B 3 (x, y, z) = 0 is a maximum for B 3 (x, y, z). 

Because the data around the extrema don’t change sign, it is a challenge 
to find these points. Fortunately, the gradient components around extrema 
does change sign. Let F'(x,y,z) be a 3-D function. Let (xo,yo,2o) be an 



extremum of F then the following conditions exist: 


dF(xo,yo,zo) _ q 
dx 

dF(x 0 ,y o 1 zo} _ q (25) 

dy v 7 

dF(x 0 y yo y zo) _ q 
dz 

Our problem has been then reduced to finding the points that satisfy one of 
the following conditions: 

VBi = 0 

VB 2 = 0 (26) 

VB 3 = 0 


Marching cube 2 is a high-resolution 3D surface construction algorithm. 
It is very helpful for us to locate zero gradient points. First, we generate 
isosurfaces of: 


si = = 0 

52 = dF{ iy' - z} = 0 

53 = = 0 


(27) 


then, we derive the intersections of these three surfaces. 

Now we find the locations where the gradients of B x ( x,y,z ), B 2 (x,y,z) 
and B 3 (x, y, z) are zero, but it is still possible that these places are just a local 
maximum or minimum for B \ , B 2 and i? 3 . We check them by substituting 
these values back into B x , B 2 and B 3 . If these functions are actually zero, 
then they are the real points. 


Application And Future Work 

Application 

As an application, we studied the Boussinesq problems with 2 forces. The 
results are quite encouraging; we found points of both double and triple 
degeneracy. All points of triple degeneracy are isolated from one another; 
however, the lines and surfaces formed from the points of double degeneracy 
extend to these isolated points of triple degeneracy. Although we studied 
only one case, our findings point to the possibility that points of double 
degeneracy in 3-D space can’t be isolated. Since the points of triple degen- 
eracy are what we are really interested in, we then analyze the topologi- 
cal structure at a small neighborhood surrounding a degenerate point where 
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\ x = A 2 = A 3 . We use hyperstreamlines to represent our eigenvector fields. In 

the figures shown below, forces are acting along the ^-coordinate direction. 
(See Plate 1.) The tensor field surrounding the degenerate point is repre- 
sented by hyperstreamlines. 1 The hyperstreamlines are integrated along the 
major eigenvector where the cross section is comprised of the two transverse 
eigenvalues. A typical cross section is an ellipse with the major and minor 
axes scaled by the transverse eigenvalues. When the eigenvalues of the corre- 
sponding transverse eigenvectors are not zero, the cross section is an ellipse 
and thus the hyperstreamlines appear as a tube; when one of the eigenvalues 
is zero, the cross section reduces to a line and if it remains zero along the 
trajectory the hyperstreamlines will appear as a ribbon; when both of the 
eigenvalues are equal, the cross section becomes a circle (this is a particular 
double degeneracy case) and if they remain equal along the trajectory then 
the hyperstreamline follows a line of double degeneracy. The color of the 
hyperstreamlines encode the magnitude of the major eigenvalue. If we are 
very close to the degenerate point, the value changes dramatically during a 
very short distance, therefore we linearize the log of the eigenvalues in our 
color mapping (except for Plate 1, which is taken from previous work). 

Plates 1 through 4 present the preliminary results from a Boussinesq 
problem with two forces. 

Plate 1 is a global view of hyperstreamlines of the major eigenvector field. 
Frame box dimensions range from -2.0 to 2.0 in each direction. Forces (repre- 
sented by two red arrows) are acting at (—1.0, 0.5, 0.0) and (—1.0, —0.5, 0.0) 
and are pointing toward +x direction. The two points of triple degener- 
acy (represented by two red dots) are at (1.05,0.0,0.5) and (1.05,0.0, 0.5). 

Color encodes the magnitude of the major eigenvalues. 

Plates 2 through 4 are three different points of view of an enlarged area 
around the degenerate point (1.05,0.0,0.5) in Plate 1, where x ranges from 
1.042 to 1.062, y ranges from -0.008 to 0.012 and z ranges from 0.492 to 
0.512. 

From Plate 2, we can see that the cross section of a hyperstreamline along 
the x direction (the same direction as the force) is a circle which means the 
two transverse eigenvalues are equal. 

From Plates 2 and 3, we can see that a hyperstreamline in the x direction 
propagates along the force first, then turns right before the triple degenerate 
point. The hyperstreamlines behind the degenerate point propagate in the 
planes perpendicular to the force where one of the other two eigenvalues is 



The pattern of the hyperstreamlines shown above represents only one type 
of triple degeneracy. The major eigenvectors propagate along the forces until 
they are very close to the degenerate point, then they all stay in the planes 
perpendicular to the forces. The pattern of their propagation in these planes 
is very similar to hyperstreamlines in a 2-D tensor field with a trisector point. 
The separatrix with the same direction as two forces is also a line comprised 
of points of double degeneracy for the other two eigenvalues, therefore the 
cross section of the hyperstreamline along this separatrix is a circle. The 
hyperstreamlines that propagate in the planes perpendicular to the forces 
are more like ribbons because one of the medium or minor eigenvalues are 
zeros in those areas. 

Future Work 

The results presented in the previous section reveal only a small part of the 
complex nature of 3-D tensor fields. The research planned for the coming 
year, based on the preliminary findings and in accordance with the original 
proposal, will focus on: 

1. Extension of topological rules from 2-D to 3-D tensor fields. 

2. Identification and classification of various types of points of double and 
triple degeneracy. 

3. Developing technique to locate continuous lines of double degeneracy. 

4. Studying the various topological features, i.e. lines and surfaces of dou- 
ble degeneracy, points of triple degeneracy, and the relations between 
them. 
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